Static NLO susceptibilities: 
testing approximation schemes against exact results 
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The reliability of the approximations commonly adopted 
in the calculation of static optical (hyper)polarizabilities is 
tested against exact results obtained for an interesting toy- 
model. The model accounts for the principal features of typi- 
cal nonlinear organic materials with mobile electrons strongly 
coupled to molecular vibrations. The approximations intro- 
duced in sum over states and finite field schemes are ana- 
lyzed in detail. Both the Born-Oppenheimer and the clamped 
nucleus approximations turn out to be safe for molecules, 
whereas for donor- accept or charge transfer complexes devi- 
ations from adiabaticity are expected. In the regime of low 
vibrational frequency, static susceptibilities are strongly dom- 
inated by the successive derivatives of the potential energy 
and large vibrational contributions to hyperpolarizabilities 
are found. In this regime anharmonic corrections to hyperpo- 
larizabilities are very large, and the harmonic approximation, 
exact for the hnear polarizability, turns out totally inadequate 
for nonlinear responses. With increasing phonon frequency 
the role of vibrations smoothly decreases, until, in the an- 
tiadiabatic (infinite vibrational frequency) regime, vibrations 
do not contribute anymore to static susceptibilities, and the 
purely electronic responses are regained. 



I. INTRODUCTION 

Photonic and optoelectronic applications heavily rely 
on materials with high nonlinear optical responses, and 
the development of new and more efficient materials for 
nonlinear optics (NLO) is a key issue in present days re- 
search. In this context, organic materials play an im- 
portant role in view of their large and fast NLO re- 
sponses and ofpthe tunability of their properties via chem- 
ical synthesisJjQ In order to trace guide-lines to or- 
ganic synthesis, much effort is devoted to understand 
the origin of nonlinearity in these systems and to re- 
late NLO pcesponses to electronic structure and molecular 
geometry.^ Q Mobile electrons are required for large non- 
linearity and conjugated molecules and/oij-,pplymers are 
materials of choice for NLO applications.0'cl|-yibrations 
are strongly coupled to conjugated electronsiB the effect 
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of electron-phonon (e-ph) coupling in NLO resncmses at 
optical frequencies is not fully understood yeiDEj, but 
general agreement has emerged on the, importance of e- 
ph coupling in static NLO responses .t5il3 

Static susceptibilities are the opposite of the successive 
derivatives of the ground state (gs) energy with respect 
to a static electric field (F): 



£ = £o- fiF 
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(1) 



where £o is the gs energy at zero field, ^ is the perma- 
nent dipole moment, a is the linear polarizability, and 
(3 and 7 are the first and second hyperpolarizabilities, 
respectively. As it is always the case for energy deriva- 
tives, two diffcimit methods can be adopted to evaluate 
susceptibilities .lia The first method, the so-called finite- 
field (FF) approach, basically relies on direct numeri- 
cal derivation techniques. The relevant system is the 
molecule in the presence of an external static electric 
field .£3 The corresponding Hamiltonian is diagonalized 
to calculate the gs energy (or dipole moment) at dif- 
ferent F values, and (hyper)polarizabilities are calcu- 
lated as numerical derivatives of the energy (or dipole 
moment) with respect to the electric field. The second 
approach treats -F as a perturbation on the molecular 
states, leading to the well-known sum over states (SOS) 
formulas.ll3 The main advantage of SOS over FF is that 
only unperturbed {F = 0) eigenstates are required. The 
main disadvantage is that all eigenstates and transition 
dipole moments are required, whereas FF only requires 
the gs energy (or dipole moment). Of course the two 
approaches are perfectly equivalent if exact energies and 
eigenstates, including vibrational contributions, are in- 
serted into equations. But for molecules and, more gen- 
erally, for molecular materials, exact vibronic eigenstates 
are hardly accessible, and several approximations are in- 
troduced. Approximated calculations based on FF or 
SOS approaches are clearly different and it is important 
to investigate the merit and reliability of the various ap- 
proximation schemes. 

The first approximation usually introduced is the adi- 
abatic (Born-Oppennheimer, BO) approximation. This 
is an almost ubiquitous approximation in molecular and 
condensed matter physics: accounting for the differ- 
ent masses of electrons and nuclei (and then for the 
different time scales of their motions), BO approxima- 
tion separates electronic and vibrational degrees of free- 
dom, so that each state is factorized into a product 
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of an electronic and a vibrational wavefunction. This 
approximation is fairly safe: deviations from the adia- 
batic regime are expected only for degenepie or quasi- 
degenerate stateSj-Jor narrow-band solids ,E3 or near to 
phase transitions 

In SOS approaches BO factorization is used to sepa- 
rate electronic and vibrational contributions to suscep- 
tibilities. Electronic contributions are those where all 
summations run on vibronic states in the excited elec- 
tronic manifolds, whereas vibrational contributions have 
at least one of the indices in the summation runmag over 
excited vibrational states within the gs manifold.E3 In the 
same spirit as BO, clamped nucleus (CN) approximation 
has been proposed to calculate SOS pcpressions for sus- 
ceptibilities in a two-step procedure.Ej In the first step 
the electronic (hyper)polarizabilities are calculated with 
nuclei clamped at their equilibrium positions. The re- 
sulting electronic contributions to (hyper) polarizabilities 
only contain sums over purely electronic states. In the 
second step the nuclear motion on the gs potential en- 
ergy surface (PES) is accounted for, and vibrational con- 
tributions to susceptibilities are obtained. The adiabatic 
approximation, as directly applied to SOS or within the 
CN formalism, makes the calculation of susceptibilities 
possible, provided that electronic and vibrational eigen- 
states are known. However, in view of the anharmonicity 
of the relevant PES, calculating the vibrational states 
in all electronic manifolds (as required in BO-SOS for- 
malism) or just in the gs manifold (as required in CN- 
SOS) is not easy. Often SOS calculations of static (hy- 
per)polarizabilities of molecules and/or pob£rp£r|S_are car- 
ried out in the harmonic approximationjlj'Ej lj where 
the solution of the vibrational problem can be written in 
closed form. Harmonic schemes work in the hypothesis 
that the electronic system responds linearly to the nu- 
clear motion: harmonic approximation is bound to fail 
in systems where vibrations are coupled to mobile elec- 
trons, i.e. . t|0 electrons characterized by large nonlinear 
responses. E§E3 In fact large anharmonic corrections to 
NLO responses have been recognized in recent ab wipQ 
calculations on conjugated molecules and polymersOEa 

When applied to FF approaches, BO quite naturally 
leads to a two-step procedure. First of all the adiabatic 
electronic Hamiltonian (i.e. the total Hamiltonian minus 
the nuclear kinetic energy, KE, term) is solved to get the 
gs electronic energy, that represents the potential energy 
(PE) for the nuclear motion. Its derivatives on F measure 
the PE contributions to susceptibilities. The calculation 
of the PE contribution is easy since only the gs electronic 
energy is required. Moreover, if nuclei are allowed to re- 
lax in response to F, susceptibilities calculated as PE 
derivatives naturally include vibrational contributions, 
and fully account for anharmonicity. The second step ac- 
counts for KE contribution: this is a more difficult step 
because the energy of the lowest vibrational state has to 
be calculated as a function of F, a non-trivial task due to 
anharmonicity. Moreover, in systems with large e-ph cou- 
pling, like conjugated polymers and/or oligomers, it can 



be difficult to get reliable calculations of F-dependent ge- 
ometries, and numerical instabilities of the energy deriva- 
tives are sometimes observed. E3 

In this paper we critically review the different approx- 
imation schemes for the calculation of static NLO re- 
sponses, when applied to molecules of interest for NLO 
applications. In these molecules, in fact, electron-phonon 
coupling has particularly large and non-trivial effects, 
shedding doubts on the reliability of approximations that 
work well in more conventional systems. In particular we 
will focus on an interesting toy-model for mobile electrons 
coupled to vibrations: the Holstein DA dimer. The so- 
called DA dimer is a two-state model originally proposed 
by MuUiken to describe electron donor-acceptor charge 
transfer (CT) complexes in soltttion.Ey Later on it was 
applied by Oudar and ChemlaEil to describe push-pull 
chromophores where the donor and acceptor moieties are 
linked by a 7r-conjugatcd bridge to form a molecule. The 
extension of the model to include Holstein caupling to vi- 
brations was first suggested by Person et alE3 to describe 
vibrational spectra of CT complexes, and was extensively 
discussed in the 80 's to describe many properties (rang- 
ing from viiajatipnal spectra to structural instabilities) 
of CT saltsrJL^ Holstein coupling in the DA dimer has 
also been adopted to describe vibrations amd/or polar sol- 
vation effects in push-pull chromophorealj: an extensive 
comparison with spectral data for several key molecules 
proved that this model contains the mai n ingre dients to 
understand nonlinearity in these systems .EjEjO In spite 
of its richness, the model is very simple and in the next 
Section we will describe its exact, non-adiabatic eigen- 
states, as well as its BO solution. Sections HI and IV 
discuss the approximations usually introduced in SOS 
and FF approaches, respectively; Section V summarizes 
main results. 



II. THE HOLSTEIN DA DIMER: 
NON-ADIABATIC AND ADIABATIC 
EIGENSTATES 

The DA dimer describes two electronic states, \DA) 
and \D'^A~), linearly coupled to a harmonic phonon, as 
depicted in Fig. 1. The relevant Hamiltonian is {Ti — 1): 

n = 2zop - V2ta^ + i (cj2g2 + p2^ _ y/2^pUjQp - fiF. 

(2) 

The first two terms in the above equation describe the 
electronic Hamiltonian, where p — {1 — a z)/ '2 is the ion- 
icity operator, and a^, represent the Pauli spin oper- 
ators. The third and fourth terms describe a harmonic 
vibration of frequency lo {Q and P are the correspond- 
ing coordinate and momentum), coupled to the electronic 
system. The strength of the coupling is measured by egp, 
the relaxation energy of \D^A~) (see Fig. la). The last 
term accounts for a static electric field, F, interacting 
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with the dipole moment operator, defined, by /t = /ioPj 
with /uo = {D+A~\fi\D+A^)M In the following, ener- 
gies are expressed in \/2t units. We underline that the 
single-mode Hamiltonian in the above equation captures 
the essential physics of e-ph coupling, at least if one is 
not intcji:cstfjd,in the detailed description of vibrational 
spcctra.lla'La'L3 The extension to the multimode case, pos- 
sible in principle, is very demanding if the non-adiabatic 
solution is required. 

The above Hamiltonian can be numerically diagonal- 
ized on the basis of the direct product of the two elec- 
tronic states \DA) and \D~^A~) and of the reference vi- 
brational states (i.e. the eigenstatesr^af the harmonic os- 
cillator in the third term of Eq. The basis is trun- 
cated by fixing a maximum number of phonon states, M; 
the corresponding 2M x 2AI matrix can be diagonalized 
up to fairly large M values, yielding numerically exact 
non-adiabatic eigenstates. The minimum M required to 
get convergence depends on the model parameters and 
on the properties of interest. All results presented in the 
following have been obtained with M = 50. 

The exact non-adiabatic eigenstates at F = 
are inserted into standard SOS expressions for static 
susceptibilitie£3 to get the exact NLO responses re- 
ported as continuous lines in Fig. 2 for egp = 1 and 
a few Lo values. By the way, exactly the same results 
are obtained, within FF approach, from the successive 
derivatives of the exact non-adiabatic gs energy vs the 
applied field. In the same figure, dot-dashed lines cor- 
respond to the bare electronic susceptibilities, Xo"^ 
to the response of the two-state model with no e-ph cou- 
pling. The deviations of the continuous lines from the 
dot-dashed lines measure the vibrational contribution to 
static responses. Dashed lines show the PE contribution 
to susceptibilities. Analytical expressions for these curves 
have already been reported in RefJlS: PE-susceptibilities 
are independent of uj. Fig. 2 clearly shows the evolu- 
tion of static susceptibilities with phonon frequency: in 
the low-w limit, KE contributions vanish, and the exact 
curves tend to the limiting PE results. Vibrational con- 
tributions to static susceptibilities are very large in this 
limit, and increase with the order of nonlinearity. With 
increasing uj, the vibrational contributions to static NLO 
responses decrease: in the antiadiabatic limit {uj oo) 
the vibrational contributions to static hyperpolarizabili- 
ties vanishes, and the exact curves tend to the bare elec- 
tronic responses. 

For real molecules or complexes, the large number of 
electronic states makes non-adiabatic calculations very 
demanding, and BO approximation is usually invoked. 
In BO approximation the effective electronic Hamilto- 
nian, Hei = TC — 1/2P^, is defined by subtracting the 
nuclear KE term from the total Hamiltonian. In our 
case Hei describes two electronic states separated by an 
energy gap that linearly idepends on Q, according to 
2z{Q) = 2zq — yj/2espUjQ&B By diagonalizing the elec- 
tronic 2x2 matrix one gets analytical expressions for 



the Q-dependent ground and excited state energies, as 
shown in Fig. lb. It is important to recognize that, even 
if the Hamiltonian in Eq. ^ is defined in terms of har- 
monic reference states (Fig. la), as a consequence of e-ph 
coupling, the ground and excited state PES are anhar- 
monic (Fig. lb). The anharnionicity of the potential pre- 
vents the analytical solution of the vibrational problem 
on either the ground or excited state; however numeri- 
cally exact vibrational states can be calculated in both 
manifolds. In particular, the eigenstates of the harmonic 
oscillator with frequency to, centered at the relevant equi- 
librium position, are a good basis for the vibrational 
problem on either PES. The corresponding vibrational 
Hamiltonian is the sum of a KE term, whose matrix el- 
ements are trivial in the adopted basis, plus a PE term, 
whose matrix elements are calculated via numerical inte- 
gration. Of course the vibrational matrix is diagonalized 
on a basis truncated to a large enough number of phonon 
states as to get convergence. Results presented in this pa- 
per have typically been obtained with 20 phonon states. 
Once BO eigenstates are obtained, the (transition) dipole 
moments entering SOS expressions can be calculated via 
numerical integration. 

Static susceptibilities calculated within BO approxi- 
mation are indistinguishable (in the scale of Fig. 2) from 
the exact ones, as long as u; < 0.2. Fig. 3 compares non- 
adiabatic (continuous lines) and BO (dashed lines) static 
susceptibilities for uj = 0.5, where deviations appear. Of 
course BO approximation becomes worst with increas- 
ing UJ and is totally untenable for uj > 1. The adopted 
value of the small polaron binding eneiay. e_, ^ \/2t ap- 
plies to both push-p|tt]|LchromophoresE3LEl~L3 and to CT 
complexes and salts.a'EJ In push-pull chromophores the 
typical value of \f2t ~ 1 eV is much largejp-tban typical 
vibrational frequencies [uj ^ 1000 cm~^)E3E3"E3 and the 
dimensionless uj ^ 0.1 — 0.2 is safely in the BO regime. 
In CT salts, instead, ^/2t is much smaller (- 0.2 eV)^^^ 
and the corresponding dimensionless w ~ 0.5 — 1 suggests 
seizable non-adiabatic effects in NLO responses of these 
(narrow-band) systems. 



III. SOS SUSCEPTIBILITIES 

The separation of electronic and vibrational degrees 
of freedom in BO approximation not only leads to sim- 
pler calculations, but also allows for additional insight. 
When applied to SOS expressions, BO quite immediately 
leads to a separation of electronic and vibrational contri- 
butions. If \Rv) indicates a BO vibronic state, i.e. the 
product between an electronic state and a vibrational 
state \v) in the R majufold, the SOS expressions for static 
susceptibilities readO 



(OG|/i|i?w)(wi?|/i|GO) 
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(OG\fi\Rv){vR\Ti\Su){uS\fi\GO) 



R,v,S,u 
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where the gs, |G0), is excluded from the summations, 
cjj^y is the frequency of the |G0) \Rv) transition, and 
Jl — fi — {OG\fi\GO). In the two-state model, the sum- 
mations on R, S and T only run on the two electronic 
states |Gl-|and \E). Then, according to the standard 
definitioncEl, the electronic contribution to susceptibili- 
ties is given by terms in Eq. with R = S = T = 
and the vibrational contribution is described by terms 
where at least one of R, S or T states coincides with |G). 
We underline that electronic contributions to BO-SOS 
susceptibilities do not coincide with the responses of the 
bare electronic system (dot-dashed lines in Fig. 2). In 
fact summations in Eq. run on vibronic states and 
not on pure electronic states. In the reliability range of 
BO approximation, closure on vibrational states in the 
excited state manifold works fine in reducing sums on in- 
termediate vibronic states into sums on pure electronic 
states, but the property is evaluated for the true vibronic 
gs (I GO)) and not as a pure electronic property (mediated 
over |G)): the electronic contribution to BO-SOS sus- 
ceptibilities partly include zero-point vibrational average 
(ZPVA) corrections, as defiiipdia Ref.Ell. 

The CN approximationoo'EJ introduces a different 
way of partitioning electronic and vibrational degrees of 
freedom. In this approach the electronic contribution to 
the n-th order susceptibility is calculated as the n-th or- 

(n) 

der bare electron susceptibility, ;\;q . Two more contribu- 
tions are then added to get the total response: (a) ZPVA 
corrections, that account for the difference between the 
electronic susceptibility calculated at the bottom of the 
gs PES (i.e. at the equilibrium geometry) and that rele- 
vant to the ground vibronic state (due to anharmonicity 
the relevant geometry does not coincide with the equilib- 
rium one). These corrections are usually calculated via 
an expansion of electronic susceptibilities on (5;lij (b) vi- 
brational contributions, that are calculated in terms of 
summations running on the vibrational states in the gs 
manifold only, as detailed in Ref.cil. 

BO-SOS and CN-SOS are based on the same adiabatic 
approximation and lead to basically identical results at 
least as long as the adiabatic approximation itself is re- 
liable (approximately w < 0.5). However, the partition- 
ing of the susceptibilities is different and there is not a 
one-to-one correspondence between terms appearing in 
the two schemes: CN-SOS kills a few terms appearing 
in BO-SOS, and, at the same time, some of tJae terms, 
that are usually considered as corresponding,^ do not 



have exactly the same meaning in the two approaches. 
As a matter of fact, differences between CN and BO re- 
sults amount to ^ 7% and ^ 15% in the electronic /? and 
7, respectively, and ~ 2% and ^ 3% in the vibrational 
contributions to /3 and 7 (these estimates refer to the pa- 
rameters Esp = I, Lo — 0.2; for lo = 0.1 deviations are 
even smaller). These small differences compensate each 
other, leading to negligible (< 1%) differences in the total 
(electronic-l- vibrational) responses. 

Our numerical solution of the adiabatic problem for 
the two-state model allowed us to test BO and CN ap- 
proximations without introducing additional approxima- 
tions. However getting exact vibrational eigenstates is 
hardly possible for complex (molecular) structures. For 
this reason CN-SOS calculations are often implemented 
by invoking theJjiarmQiiic approximation for the vibra- 
tional problem. 0000 As already stressed in the In- 
troduction, the harmonic approximation fails in systems 
where vibrations are coupled to electrons with large non- 
linear responses. The relation between anharmonicity 
and nonlinearity is very well apparent in the adopted 
model, where both Q and F variables are coupled to 
the same electronic operator, p (cf Eq. (^). Successive 
Q-derivatives of the ground and excited state PE {£c 
and £e, respectively) are therefore directly related to F- 
derivatives, i.e. to the bare electronic susceptibilities, 
Xo"\ as follows: 



G/E 



9Q" 



1/2, ,nXo 



-1) 



Mo 



(4) 



where S2n is the Kronecker-(5, equal to 1 if 71 = 2, zero oth- 
erwise. From this equation it turns out that systems with 
large /ij/perpolarizabilities, i.e. with large Xo^\ for n >2, 
are characterized by largely anharmonic ground aruLcx- 
cited state PES. As it has been discussed in RefHjO, 
the anharmonicity of vibrations coupled to delocalized 
electrons is hardly detected in coherent spectral measure- 
ments (e.g. electronic and vibrational absorption and/or 
Raman spectra), where experimental data can be sat- 
isfactorily reproduced in terms of parabolic PES with 
effective curvatures, but it shows up with large effects 
in incoherent spectral measurements (e.g. steady-state 
and/or time-resolved emission) £3 Large anharmonic cor- 
rections cii£_also expected in static nonlinear optical 
responses, as it is confirmed., by extensive ab initio 
results on conjugated materials .EZl'ta 

The dotted lines in Fig. 2 show the (cj-independent) 
susceptibilities |Galculated in the best harmonic approx- 
imation (BHA)E3. In this approach the ground and ex- 
cited state PES are approximated by the two parabolas 
that best fit the two anharmonic PES at the gs equilib- 
rium position (i.e. at the minimum of the gs PES): the 
relevant curvatures are then defined by Eq. (^ in terms 
of the bare electronic linear polarizability evaluated at 
the equilibrium geometry. Moreover, all electronic prop- 
erties appearing in CN-SOS expressions are truncated 
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to the linear term in Q.c2l We defer a detailed compar- 
ison of BHA and exact curves to the next section, here 
we only underline that, as far as a is concerned, BHA 
and PE estimate exactly coincide, confirming that in the 
to regime anharmonic corrections to the linear polar- 
izability tend to vanishE3'Lj. In the same limit, nonlinear 
susceptibilities (/? and 7) are largely amplified by the an- 
harmonicity of the gs PES, as demonstrated in panels (a) 
and (b) of Fig. 2, by the large deviations of BHA results 
(dotted lines) from PE curves (dashed lines). 

IV. FINITE-FIELD SUSCEPTIBILITIES: 
POTENTIAL AND KINETIC ENERGY 
CONTRIBUTIONS 

In SOS approaches, susceptibilities are naturally par- 
titioned into electronic and vibrational contributions. 
Within FF approach, on the opposite, PE contributions 
to the susceptibilities are naturally separated from con- 
tributions due to the nuclear KE. As already discussed in 
the Introduction, the lowest eigenstate of the electronic 
Hamiltonian as defined in the adiabatic approximation, 
is the PE for the motion of nuclei (in the gs manifold, 
of course) . The total gs energy is obtained by summing 
the nuclear KE to the PE. Then, as long as BO applies, 
susceptibilities, i.e. the successive derivatives of the gs 
energy with respect to an applied electric field, can be 
calculated as sums of PE and KE _F-derivativcs. As a 
matter of fact, one can further separate the PE contribu- 
tion to susceptibilities into an electronic and a vibrational 
part by simply comparing derivatives calculated at fixed 
nuclear geometry with those taken by allowing nuclei to 
relax following the application of the electric field. Ej This 
separation is however fairly artificial and does not add so 
much to our understanding. 

The nuclear KE vanishes in the lu = limit, and the 
(w-independent) PE susceptibilities, reported as dashed 
lines in Fig. 2, represent the zero- frequency limit of the 
exact susceptibilities. Analytical expressions for the PE- 
susceptibilities p£ the Holstein DA dimer were already 
reported in Ref.EJ. More generally, the calculation of PE- 
susceptibilities is easily implemented in quantum chem- 
istry calculations, in fact it only requires the gs energy 
calculated at the relaxed geometry for different values 
of an externally applied field. Of course, geometry op- 
timization is a crucial step in FF approach to hyperpo- 
larizabilities and it is extremely important to fully re- 
lax the geometry for the chosen molecule in order to 
avoid spurious 'straial or 'instability' contributions to 
hyperpolarizabilities .E3 

The calculation of KE contributions is more difficult, 
since the F-dependence of the lowest vibrational state 
in the anharmonic gs PES is needed. In the adopted 
model, however, the calculation is feasible and, as long 
as BO applies, leads to basically exact results. Specif- 
ically, much as it occurs for BO-SOS and CN-SOS, for 



LU < 0.2 exact (continuous) curves in Fig. 2 are indis- 
tinguishable from BO-FF results. The nuclear KE con- 
tributes to susceptibilities in two different ways. First 
of all, due to anharmonicity, the molecular geometry in 
the vibronic gs is different from the equilibrium geome- 
try (corresponding to the minimum of the gs PES). The 
PE susceptibilities have, in principle, to be corrected to 
account for this effect. The correction is however very 
small, and we found it always negligible in the investi- 
gated parameter range. The second contribution stems 
from the _F-dependence of the nuclear KE itself: it is this 
contribution that indeed accounts for the deviations of 
the exact curves form the (dashed) PE curve (at least in 
the BO regime, u; < 0.2, where non-adiabatic corrections 
are negligible) . KE contributions are of course very small 
for low uj, but they increase with increasing lu, leading 
to an overall decrease of the vibrational amplification of 
the static NLO responses. This is by no means acciden- 
tal: in the antiadiabatic limit (cu ^ V^t) phonons cannot 
contribute to static susceptibilities and, with increasing 
LU, KE contributions progressively increase to counterbal- 
ance PE contributions. 

KE contributions to susceptibilities exactly vanish in 
the harmonic approximation. In fact, the equilibrium 
position in any harmonic vibrational state coincides with 
the bottom of the PES; moreover the nuclear KE is pro- 
portional to the harmonic frequency, i.e. to the curvature 
of the PES. For a parabolic PES, this quantity is obvi- 
ously independent of Q, and hence of F (cf Eq. (|^)). We 
are then in the position to prove that the harmonic ap- 
proximation is totally uncontrolled for the calculation of 
/lypeTpolarizabilities in any frequency regime. At low lu, 
where PE contributions dominate, the harmonic approx- 
imation cannot account for anharmonic PES: this has no 
effect on the linear response, and the BHA estimate, of a 
coincides with the exact curve in the uj — limit .E3 An- 
harmonic contributions to PE-ZiT/perpolarizabilities are 
instead large, as demonstrated by the large deviations 
of BHA curves from the lu ~ limiting curves. At the 
same time, not accounting for KE corrections, harmonic 
results do not represent a reliable approximation in the 
large lu regime either. It is possible that in a given pa- 
rameter range BHA gives good estimates of a specific 
response, but this is due to an accidental compensation 
of two large errors. The harmonic approximation is un- 
reliable in the calculation of /lyperpolarizabilities, and, as 
it turns out clearly from Fig. 2, it is not possible to devise 
any parameter range where harmonic results represent a 
good approximation for all static susceptibilities. 

V. CONCLUSIONS 

The MuUiken DA dimer model, extended to account 
for the Holstein coupling to vibrations (or equivalently to 
other slow degrees of freedom, including the orientational 
solvation coordinate), contains the essential physics to 
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understand the subtle interplay between electrons, vi- 
brations and external perturbations (namely an electijicj 
field) governing NLO responses of molecular materials.E-3 
This simple toy-model is then very useful to test the 
reliability of several approximation schemes usually in- 
troduced in more refined quantum chemical description 
of these materials. It is important to recognize that 
some very well-known and very widely applied approx- 
imations can fail when applied to model the properties 
of materials with strongly nonlinear behavior. Pertur- 
bation theory, and particularly linear perturbation the- 
ory, is clearly inadequate in this context, and we already 
proved its inadequacy in reprO|d«eiM spectral proper- 
ties of push-pull chronnophoresEaEaO as well as their 
static susceptibilities.t3 More subtle effects have been 
demonstrated in two-photon absorption (TPA) spectra, 
where the standard Condon approximation badly fails 
even when it properly reproduces one-photon absorption 
spectra.Cj Even more interestingly, an important vibra- 
tional channel is found to contribute to TPA spectra, in 
addition to standard electrcauc channels, leading to large 
effects in observed spectraJlJ The harmonic approxima- 
tions works fine in TPA spectra and more generally in 
vertical (coherent) processepyibut it is inadequate to de- 
scribe incoherent processes.E3 

In this paper we have addressed the approximations 
usually introduced in the calculation of static NLO re- 
sponses. The adiabatic approximation is fairly safe, at 
least as long as phonon frequencies are not too high (ap- 
proximately uj < 0.5 1, cf Fig. 3). It can be applied in 
different ways, BO or CN, in SOS approaches leading to 
a slightly different partition of electronic and vibrational 
contributions to susceptibilities, but to the same overall 
result. In FF it leads to a partitioning of potential and 
kinetic energy contributions to susceptibilities, but oth- 
erwise leads to the same result as in SOS. Closure over 
vibrational states is also a very good approximation, as 
long as the adiabatic approximation applies. 

Potential energy contributions are clearly domi- 
nant over kinetic energy contributions in the low- 
vibrational frequency regime. The anharmonicity of 
the gs PES therefore gives very large corrections to hy- 
perpolarizabilities: the harmonic approximation is unten- 
able in this regime. With increasing phonon frequency, 
kinetic energy contributions globally reduce vibrational 
corrections until, in the antiadiabatic limit (u — > oo) only 
electrons respond to static fields to regain purely elec- 
tronic susceptibilities. Kinetic energy corrections vanish 
in the harmonic approximation and this approximation 
fails in the high frequency regime too. 

Both FF and SOS approaches can be easily imple- 
mented in quantum chemistry calculations: the former 
has the advantage of only requiring knowledge about the 
lowest eigenstate, whereas all eigenstates are required in 
SOS. However only F = eigenstates are needed in SOS, 
whereas the -F-dependence of the gs energy (or dipole mo- 
ment) is required in FF. Anharmonicity is important and 
must be accounted for: in the low-vibrational frequency 



limit, where kinetic energy contributions are negligible, 
FF calculation is very convenient and fully accounts for 
anharmonicity. In the intermediate frequency regime the 
PE-FF calculation overestimates the vibrational amplifi- 
cation of NLO responses: either KE corrections have to 
be introduced in the FF scheme, or one has to make re- 
sort to SOS calculations, but in any case anharmonicity 
has to be accounted in order to get reliable estimates of 
/i?/pe?polarizabilities. 
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20 
Q 

FIG. 1. Potential energy surfaces for (a) the basis states {\/2t = 0) and (b) the exact eigenstates {\/2t = 1) of the Hamiltonian 
in Eq. (1). All curves are calculated for zo = 1.3 and tsp = 1. 
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FIG. 2. Static (hyper)polarizabilities as a function of the 
ground state ionicity pgs, for e^p = 1 and different uj val- 
ues (\/2t units). Dashed hnes: potential energy contribution 
to susceptibilities, corresponding to a; = limit (see text); 
dot-dashed lines: bare electronic susceptibilities, correspond- 
ing to the UJ ^ oo limit (see text); continuous lines report ex- 
act susceptibilities calculated for ui = 0.05, 0.2, 1.0, smoothly 
evolving from the w = to the oi — » oo limits; dotted lines 
report susceptibilities calculated in the best harmonic approx- 
imation (see text). In panel (a), dotted and dashed lines are 
exactly superimposed. For the calculation of a, (3 and 7, 
dipole moments have been expressed in //o units. 




FIG. 3. Static (hyper)polarizabilities as a function of the 
ground state ionicity Pgs, for tsp = 1 and w = 0.5 {V2t 
units). Continuous lines show non-adiabatic results; dashed 
lines show BO results. 
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